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We present a functional renormalization group flow for many-fermion lattice models into phases 
with broken spin-rotational symmetry. The flow is expressed purely in terms of fermionic vertex 
functions. The symmetry breaking is seeded by a small initial anomalous self-energy which grows 
at the transition scale and which prevents a runaway-flow at nonzero scales. Focusing on the case 
of commensurate antiferromagnetism we discuss how the interaction vertex can be parametrized 
efficiently. For reduced models with long-range bare interactions we show the results of standard 
mean-field theory are reproduced by the fRG and how anisotropics in the spin sector change the 
flows. We then describe a more efficient decomposition of the interaction vertex that should allow 
for the treatment of more general models. 
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I. INTRODUCTION 

Spontaneous symmetry breaking is an important 
building block of our understanding of interacting mat- 
ter. Typically, the symmetry breaking occurs at energy 
scales much smaller than those at which the correspond- 
ing theoretical models are formulated. This makes a 
quantitative theoretical description of these models and 
hence of the respective materials challenging. In a func- 
tion renormalization group (fRG) approach (for a recent 
review see Ref. Q]), flow equations interpolate between 
the parameters of the microscopic (bare) action at an 
high energy scale and those of the one-particle irreducible 
functional, which can be regarded as an effective action 
at lower scales. 

For the description of spontaneous symmetry break- 
ing in correlated electronic systems, one often decouples 
the two-particle interaction via a Hubbard-Stratonovich 
transformation. Then the fermions can be integrated 
out, and the remaining bosonic theory of the composite 
field(s) can possibly, as a first approximation, be treated 
in saddle-point approximation. However, in certain rele- 
vant cases such as unconventional superconductivity, the 
bare interaction does not contain the necessary attrac- 
tive component in the desired channels, and a more suit- 
able effective interactions needs to be constructed, e.g. 
by fermionic perturbation theory. Furthermore, even if 
this first problem is not present, the mean- field solu- 
tion ignores potentially important corrections that can 
lead to change in the resulting energy spectrum, or to 
the renormalization or even disappearance of the order 
due to collective fluctuations. Also, integrating out the 
fermions completely can lead to severe complications, e.g. 
due to diverging coupling constants of the higher order 
terms 2 . In particular for the latter aspects, it may be 
advantageous to consider the renormalization-group flow 
for a mixed system consisting of fermions and exchange 
bosons. At the onset of symmetry breaking, the renor- 
malized potential for the zero-momentum mode of such 



an exchange boson has degenerate minima correspond- 
ing to a nonzero order parameter^. In the symmetry- 
broken regime, the effective action may be expanded 
around one of these minima, giving rise to flow equations 
for the radial- and Goldstone-mode propagators and the 
fermionic gap, as in the case of a fermionic superfluid^. In 
a series of publications^—, the low-temperature physics 
of the two-dimensional (2D) Hubbard model, which may 
be regarded as a generic model for high-T c compounds, 
was analyzed within such an approach. In these pa- 
pers, the generation of an attractive ci-wave pairing in- 
teraction is taken into account by rebosonizing fermion- 
fermion interactions generated in the flow. In the par- 
tially bosonized formulation, however, the momentum 
dependence that is not captured by the bosonic channels 
has been dropped. Moreover, a Hubbard-Stratonovich 
decoupling of the bare interaction adds some bias to the 
description. For the Hubbard interaction, this problem 
has been circumvented by only bosonizing renormaliza- 
tions of the interaction^. 

While these partially bosonized flow studies have their 
merits and should be extended further, one has to pay 
some more attention in order to capture effects (such 
as the generation of attractive pairing components) that 
one gets rather straightforwardly in purely fermionic ap- 
proaches. In particular in cases with competing orders 
this might become more cumbersome. In such situa- 
tions, purely fermionic flows have been used successfully 
at scales above the spontaneous symmetry breaking (for 
a recent review see Ref. [l|). So, in another string of 
works 9 - - — one has tried to continue these purely fermionic 
flows into the symmetry breaking regime. There is, how- 
ever, a crucial difference to the mixed flows: If the bare 
action does not contain a symmetry-breaking seed field, 
the tendency toward spontaneous symmetry-breaking in 
such a purely fermionic description manifests itself in 
a flow to strong coupling, such as a Cooper or spin- 
density-wave (SDW) instability. An instability analy- 
sis of that type has — among other systems — been 
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performed for the 2D Hubbard modeli^r— and the iron 
pnictides 16 . For a non- vanishing seed field, mean- field 
models for broken translational— and U(l)2iil symme- 
try have been studied. These models are exactly solv- 
able in the thermodynamic limit, as only diagrams in 
one channel contribute. In order to reproduce the ex- 
act result, the fermionic flow needs be considered in the 
Katanin truncation- 8 . Beyond mean- field models, super- 
fluidity in the attractive 2D Hubbard model was studied 
within a purely fermionic description^ 2 -, and recently in 
channel-decomposed approach^. 

Phases of broken spin-symmetry, in contrast, have not 
been studied yet with these extended RG flows. The com- 
petition (and possible coexistence) of antiferromagnetism 
(AFM) and superconductivity (SC), however, is of prime 
interest in high-T c compounds and iron pnictides (see for 
example Refs. I20I andl2ll). Hence, fRG flows should also 
be continued into the SU(2)-broken phase. 

The goal of the present work is to lay the foundations 
for doing so in realistic models. It is organized as fol- 
lows. In Sec. HH we given an efficient parameterization 
of a spin-symmetry broken interaction and derive RG 
flow-equations for the respective coupling functions. In 
the sequel, we analyze reduced mean-field models for a 
spin-density- wave (SDW) phase. SecHUis devoted to an 
isotropic interaction. The effective action is parameter- 
ized further and couplings corresponding to the bosonic 
radial and Goldstone modes are identified. The resulting 
flow equations are found to be equivalent to resummation 
within the random-phase approximation (RPA) and the 
self-consistent gap equation. The generalization to an 
anisotropic mean-field interaction in Sec. IIVI is straight 
forward. Before we finally conclude in Sec. IVI| we pro- 
pose a channel decomposition of the flow equations in the 
SU(2)-broken phase in Sec fVl which can also be used for 
an improved treatment of more general models. 



II. GENERAL CASE 
A. Parametrization of the interaction 

In Ref. the four-point vertex was parametrized for 
theories that remain SU(2)-symmctric in spin space. The 
solution of the flow equations for appropriate cases then 
yields a flow to strong coupling where at least one class of 
components of the four-point vertex diverges at a nonzero 
critical scale. If the Fermi surface is nested, e.g. due to 
flat parallel sides, and if the initial short-range interac- 
tion is repulsive, the flow to strong coupling usually indi- 
cates an antiferromagnetic instability where the static 
susceptibility in this channel becomes infinitely large. 
The interpretation of this instability is that below this 
critical scale, antiferromagnetic order sets in. This regu- 
larizes the growth of the four-point vertex in a reasonable 
way that is discussed further below, and definitely breaks 
the SU(2)-spin rotational symmetry. Hence, in order to 
understand the low-scale regime, one needs to continue 



the flow into the broken-symmetry regime. In this section 
we address the necessary parametrization of the interac- 
tion in the SU(2)-broken phase. In order to keep the re- 
sulting numerical calculations tractable, we will restrict 
ourselves to spin-l/2-fermions exhibiting collinear anti- 
ferromagnetism, i.e. we break the spin-symmetry only in 
one direction. More specifically, we allow for (S z ) ^ 
while we still require (S x ) = (S y ) = 0. For such a 
case, a similar parametrization 2 - has been given in real 
space. The two-point Green's function only has diagonal 
entries in spin space and can then be split into a spin- 
flip symmetric and a spin-antisymmetric part G Tl ,r 2 = 
GiS Tl , T2 + G z a z Tl T2 — G Tl , with the Pauli z-matrix a z 
and the spin indices being j" or |. The elements of the 
remaining spin symmetry group are U(tp) = e %VG with 
arbitrary real if. A general two-particle interaction can 
be written as an object quartic in fermionic Grassmann 
fields, 

with generalized quantum numbers £j, containing e.g. 
fermionic Matsubara frequency, wavevector or momen- 
tum and spin quantum number along the ^-direction, Tj. 
The integration J d^i runs over all these entries. Under 
the remaining symmetry transformations U(ip), this term 
transforms according to 

implying conservation of the projected spin z- 
componcnts, t\ + r 2 = t 3 + T4. There are six spin 
configurations satisfying this constraint, namely 

tttt, uu, tin, nit, itn, mt ■ 

Since /(6, 6, 6, 6) is antisymmetric under £1 O £2 and 
6 6, the spin- dependence of the interaction can be 
parametrized using three independent momentum- and 
frequency-dependent (encoded in variables ki) coupling 
functions 

/ (6 , 6 , 6 j 6) =Vf {fa , k 2 , k 3 , k A ) &r ,ttTt 

+ Vj.(fci, k 2 , fc 3 , fc 4 ) 5t, 4,4.44. 
+ V n (k 1 ,k 2 , k 3 ,k 4 ) S T . nn 

- Vnihih, &4, h) f>i-,tut 
+ %(fc2, fei, fa, k 3 ) S T . lm 

- V n (k 2l ki,k 3 ,ki) &r,4tt4. ■ 

Due to the antisymmetry property of /(6,6,6,6), 
and V4 are antisymmetric under k± O- k 2 and k 3 O ki, 
whereas does not obey a Pauli principle constraint. 
Particle-hole symmetry implies 

Vf(kl,kj, k 3 , fei) = Vf(k4,,k3,k 2 ,kl) , 
Vj.(fcl,fe2,fe3,fel) = KUfet,fe3,fe2,fcl), 

V^{k\,k 2 ,kz,ki) = V n (k 3l k 4 ,ki,k 2 ) . 
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Restoring the SU(2) symmetry imposes the constraint 

Vf(ki, k 2 , fc 3 , ki) = V\.{k\,k 2 , k 3 , k 4 ) 

= Vxi(k 1 ,k 2 , fc 3 , ki) - V\i(k 2 , fci, fc 3 , fc 4 ) 

on the coupling functions. A global SU(2) Ward identity 
can be derived as in the U(l) case in Ref. Eq. 85. One 
obtains 

G z (ki,k 2 ) -G° z (ki,k 2 ) = ~ Jdpidp 2 dp 3 dp 4 G° z (pi,p 2 ) 
x G±(p 2 ,p 3 ) Gt(p A ,pi)V n {ki,p s ,P4,k 2 ) , (1) 



where G° z denotes the spin-antisymmetric part of the bare 
propagator. 

B. Flow equations 

This new parametrization allowing for collinear spin 
order in z-direction can now be inserted into the flow 
equations for the general two- and four-point vertices in 
Paragraph 4.1 of Ref. [22|. The simplified flow equations 
will now be given in this section. The right-hand side of 
the flow equations for the interaction can be decomposed 
into particle-particle and particle-hole diagrams 



d\V^{ki, k 2 ,k 3 , fc 4 ) = T pp (fci, k 2 ,k 3 , k 4 ) + T ph (k 4 , fc 2 , k 3 , k 4 ) . 
d\V±(ki, k 2 ,k 3 , ki) = T^ p (k 1: k 2 ,k 3 , k 4 ) + T pYl {ki,k 2 , k 3 , k 4 ) ■ 
dxVfiikx, k 2 ,k 3 , k 4 ) = T pp (fci, k 2 ,k 3 , k 4 ) + T^(ki,k 2 , k 3 , k 4 ) . 



(2) 
(3) 
(4) 



With the short-hand notation L TliT2 (p,p';q,q') — d\ [G Tl (p,p') G T2 (q,q')} for the loop diagrams, these contribution 
are as follows. 
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T t pp (fci , fc 2 , fc 3 , fc 4 ) = - dp dp' dqdq'V t (k 1 ,k 2 ,p 7 q)V t (q',p',k 3 ,k 4 )L tA (p,p';q,q'), 



(5) 



TF h (kx,k 2 ,k s ,k 4 ) = ~ dpdp' dqdq' [V t (q' ,k 2 , k 3 ,p)V t (k 1 ,p' , q,k A ) - V t (q' ,k 2 ,k 4 ,p)V t (k ly p' ,q,k 3 )} L tit {p,p';q,q') 



dp dp' dq dq' [V n (k 2 ,q',k 3l p) % (fc 4 ,p',h,q)- V n (k 2 , <?' , fc 4 , p) %(fci , p', fc 3 , q)} 
x L ui (p,p';q,q') . (6) 

Likewise, we obtain 

Tf p (ki,k 2 ,k 3 ,k 4 ) = ~ J dpdp' dqdq'Vi(ki,k 2: p,q)V i (q' : p',k 3 ,k 4 )L l ^(p,p';q : q') , (7) 

TX h (fci,fc 2 ,fc 3 ,fc4) = -J dpdp' dqdq' [V±{q' ,k 2 ,k 3 ,p)V i {k 1 ,p' ,q,k 4 ) - V_i(q' ,k 2 ,k 4 ,p)V i (k 1 ,p' ,q,k 3 )} L iA (p,p';q,q') 

-J dpdp' dqdq' [%(</, k 2 ,p, k 3 ) V n (p' , fci, g, k 4 ) - Vf±(q', k 2 ,p, k 4 ) V n (p' , ki,q, fc 3 )] 
x L t . t (p,p';q,q') 



T^(k!,k 2 ,k 3 ,k 4 ) 
T^(k!,k 2 ,k 3 ,k 4 ) 



and 

' dpdp' dqdq' V n (h,k 2 ,p,q) V n (p' , q',k 3 ,k 4 ) L t)4 .(p,p'; q,q') 

dpdp' dq dq V n (q , k 2 ,k 3 ,p) V n (ki,p',q, k 4 ) L iA {p,p ';q,q') 

dpdp' dq dq' V u {q' , k 2l p, k 4 ) Vf(fci,p', g, fc 3 ) L tA (p,p';q,q') 

^ ' "dp dp' dq dq' Vj. (<?' ,k 2 ,k 4 ,p)V n (k 1 ,p' ,k 3 ,q) L ui {p,p';q,q') 
For the flow of the self-energy, we find 

d x E t (fci , k 2 ) = - J dp dp' [S t (p, p' ) Vf (h ,p' ,p, k 2 ) - Si (p, P ' ) V n (h , p' ,k 2 , p)} 

dxEi(ki,ks) = dpdp ' [Si(p,p')V±(ki,p' \p,k 2 ) - S^(p,p')V n {p' ,h,p,k 2 )] . 



(8) 
(9) 



(10) 

(11) 
(12) 
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The first-order differential equations H H H HQ and QJ 
have to be solved together. We have not specified the 
flow parameter A yet, and these equations hold for any 
continuous parameter that appears in the quadratic part 
of the initial action only. In the following, we will use a 
sharp band energy cutoff, where the modes above scale 
A are suppressed in the free propagator. 

III. MEAN-FIELD MODEL FOR 
COMMENSURATE ANTIFERROMAGNETISM 

A. Parametrization 

Before we discuss more complicated cases, we study a 
simple mean field model for fermions on a D-dimensional 
lattice with infinitely long-ranged staggered (i.e. an- 
tifcrromagnctic) spin-spin interactions. More precisely, 
we consider lattice fermions on a .D-dimensional torus 
with circumference L in all directions and take the limit 
L — > oo. The staggering is characterized by a wave vec- 
tor Q, which does not need to be specified in full detail. 
However, for simplicity, we restrict our analysis to com- 
mensurate ordering vectors Q. Later, we will also assume 
perfect particle-hole nesting entailed by a fermionic lat- 
tice dispersion with the property — — £k+Q- This as- 
sumption is not crucial for the validity of the scheme, but 
makes the solution of the flow equations much simpler. 
The action then reads 

S = J dk rp T (k) (iko - £k) ip T (k) + — Sq • S_q , (13) 

where S Q = Jdkijj(k) <rip(k + Q) with Q = (0, Q) 
and J > 0. Note that the interaction term contains in 
total only two fc-summations, i.e. the frequency- and 
momentum-structure of the interaction is strongly or 
doubly restricted compared to the general translationally 
invariant case with three summations for momenta and 
frequencies. Its strength is renormalized by the 1 + D- 
dimensional volume ft — f3L D at temperature l/f3. In 
the following, we will consider the RG flow in the limit 
£1 — > oo which corresponds to the thermodynamic limit 
and/or zero temperature. 

In order to break both SU(2) and translational invari- 
ance, we add a spin-antisymmetric term ASq to the bare 
action. In our parametrization this corresponds to the 
following coupling functions in the ultraviolet 





2 T 

= ^[S(k 2 -k 3 - 


i-O)- 


5 (fci - k 3 -+ 


- Q)\ h^} i 




2 I 

= 1T [s(k 2 -h- 


VQ)- 


• S {kx - k 3 ■+ 


- Q)} h^} i 




2 I 

= — [25{k 2 -k 3 


+ Q)- 


+■ 6 (ki - k 3 


+ Q)\ hk.y 



Here, the superscript denotes the scale and we have in- 
troduced S(p — q) = £1 <5p , go ^p,q and 

<5{fcx} =5{k 1 +k 2 -k 3 - ki) 



as a short-hand notation for a energy and momentum 
conserving ^-function. 

Here and throughout, we use a pseudo-continuous no- 
tation, i.e. we interpret momentum and frequency in- 
tegrals Jdk as a summation fi" 1 J2k over the discrete 
momenta of the finite system at finite temperature in 
the thermodynamic limit L — > oo and send f3 — >• oo 
for zero-temperature calculations. As pointed out in 
Refs. d, [ill and [T3 particle-particle diagrams and dia- 
grams with overlapping loops vanish in the thermody- 
namic limit since the square of a ^-function produces an 
extra factor f2. In the thermodynamic limit, the only 
non- vanishing contributions to the renormalized interac- 
tion are of order f2 _1 and hence a Bethe-Salpether equa- 
tion holds 

/ A (6 , 6 , 6 , U ) = - A^ 3 ,u J dm ■■d m G x { m , V2 ) 

where the operator A antisymmetrizes the cou- 
pling function according to A^ lt ^ 2 h(^i, £2) = 
[h (£1, £2) — h (£2, £1)] /2. Note that this equation 
holds without any further assumption for the dispersion 
such as perfect nesting. Similarly, only contributions to 
the self-energy of zeroth order in f2 _1 are retained in the 
thermodynamic limit. 

For broken translational invariance, contributions to 
the interaction that violate momentum conservation are 
generated during the flow in general. In the thermody- 
namic limit, however, such interaction terms vanish due 
to the restricted momentum dependence of the bare in- 
teraction. The renormalized interaction can therefore be 
parametrized by four couplings 

V t A = J t A [6 (k 2 -k 3 + Q)-6(k 1 -k 3 + Q)] 6 {ki} /n , 

Vl = [6 (fe -k 3 + Q)-6(k 1 -k 3 + Q)] 6 {ki} /il , 

V tl = [ J *v 5 ( fc 2 - h + Q) + J X S (h - k 3 + Q)] 6 {ki} /Q . 

(14) 

In the SU(2) symmetric case, these couplings fulfill the 
constraint J-f = Jj, = J xy — J z . 

B. Flow equations 

From the ultraviolet values = J£° = 2 J, J~ = 4J 
and = 2 J these couplings flow according to 

j t = - J t 2 B t , t - J 2 z B ui , (15) 

ji = -Jf4r^,t. ( 16 ) 

Jz = ~Jz (^t.t + Ji^Ui) , (17) 
Jxy = ~J xy B^ , (18) 

with the dots denoting scale derivatives. For convenience, 
we have now suppressed the scale dependence in the no- 
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tation. The particle-hole bubble 

B Tl , T2 = ^Jdkdk' [G T1 (k,k') G T2 (k' + Q,k + Q) 
+ G Tl (k,k' + Q) G T2 {k' + Q,k)] 
corresponds to a trace in Nambu space with spinors 



*r(fc) 



ifj T {k + Q) 



(19) 



The restricted momentum dependence of the interac- 
tion leaves the momentum-conserving component of the 
self-energy unchanged during the flow. For the remaining 
anomalous 24 components, the flow equations read as 



E T (k,k') = E T 5(k - k' + Q) 
E t = -J t A t + J z Ai , 
E; = -J±A± + J z A t , 



(20) 
(21) 



with 



dk dk' S T {k,k')S(k-k' + Q) 



being the anomalous tadpole bubble. Note that if we had 
only broken the translational, but not the spin symmetry, 
the loop integrals A T and B T1>T2 in the flow equations 
Eg s . 1 1 5tf2~T1 would not depend on the spin indices and the 
number of independent couplings would be reduced to 
one according to = J^ = J xy — J z = J z . 

If the symmetry breaking term in the bare action is 
antisymmetric under a spin flip as in our case, the prop- 
erty — I] i — E-f =: A is preserved in the flow. We thus 
have Sfsf = B±,i, A\ = ~A^ and = Jj, = J z and thus 
the flow equations Eq s . [T5T [TBI and |2~T1 become redundant. 
Note that the interaction now remains invariant under 
spin flips. The flow equations to be solved finally read 



Jz = 
A = 



-2 jf-Bt-t j 
-2J z A t . 



(22) 
(23) 
(24) 



Note that the flow of the in-plane coupling J xy does not 
feed back on the flow of A and J z . This behavior is 
similar to the reduced BCS model^, where the flow of the 
amplitude vertex and the gap decouples from the flow of 
the Goldstone vertex as well. Indeed, J xy will turn out 
to account for the Goldstone modes further below. 

So far, we have not made any specific assumption for 
the dispersion. In the following, we will only consider the 
case of perfect nesting = — £k+Q in order to simplify 
the integration of the loops 25 . Since the anomalous self- 
energy E T is independent of k, we may write in Nambu 
representation 



G r (k,k') = G T (k)S(k-k') 



G T {k) 



Y 



k 2 



4 + X 2 E 2 



-iko - e k 
~XE T 



-ik + £k 



with a sharp energy-shell regulator x = ©(l e k| — A). 
The expressions for the loops are now evaluated as fol- 
lows. The Matsubara sums are performed analytically 
and Morris's lemma^ is used to evaluate Heavyside func- 
tions at their jumps. For the sake of a lean notation, we 
define 



E = \/e 2 + X 2 A 2 , 



t' = 



E 

2Tcosh ( — 
2T 



-1 -1 



where T denotes the temperature. Calculating the Mat- 
subara sum analytically, we obtain for the particle-hole 
bubble with opposite spins 



B 



d\ dkx 



t 
2E 



2E 



+ AA / dep. 

e=A,x=l "' A 



_t 

£3 ~ IP 



whereas the equal-spin bubble reads as 



B^, t ^B tA + d x dk X A 



(—-—) 

\2E 3 IE 2 ) 



A 2 a 



f tx 2 

V 2i?3 1 2E 2 J e= , =1 



-A A J dep e 



E 3 E 2 ) E 2 



TE 3 



Here, momentum integrals have been replaced by energy 
integrals with the density of states p e and ultraviolet cut- 
off W. For the tadpole bubble, Morris's lemma yields 



P\ A 



t 

2E 



If spin symmetry was restored, we would have to en- 
counter a flow to strong coupling, since then A = and 



hence B TltT2 > 0. 



In Figs. [T][3] the flow of A, J z and J xy is depicted for 
perfect nesting, zero temperature and a constant density 



of states po- At a critical scale A cr ;t oc e 



-l/(4Jp ) 



(about 



10~ 2 /po hi the case considered here), the coupling con- 
stants of the interaction grow and the gap opens. At 
lower scales, J z is then depressed by this gap while J xy 
continues to grow. The growth of J z slightly above the 
critical scale occurs as a remainder of the AFM instabil- 
ity encountered in a SU(2)-symmetric flow. If Ao is sent 
to zero, the infrared value of J z and A saturate whereas 
Jxy grows without bound in the infrared. 
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FIG. 1. Flow of the gap A for zero temperature, perfect 
nesting and a constant density of states po (W = 2/po, J = 
0.1/po). 
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FIG. 2. Flow of J z for zero temperature, perfect nesting and 
a constant density of states po (W = 2/po, J = 0.1/ po)- 
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FIG. 3. Flow of J xy according to Eqs. El El and [28] for zero 
temperature, perfect nesting and a constant density of states 
Po (W = 2/po, J — 0.1/po)- The flow equation for J xy since 
the use of Eg. 1281 avoids numerical complications. 



C. Exact solution 

The flow equations [22ll23l can be formally solved to 
reproduce the RPA result 



Jz 



2.7 



4JB tit 
4J 



1 + 4JB 



t-i 



(25) 
(26) 



The gap equation for the anomalous part of the propa- 
gator can now be derived if the tadpole is recast as 



-A B 



dkdk'G f (k,k')5{k-k' + Q) . 



Inserting Eq. [25] into Eq. [M] and integrating then yields 
the following gap equation 



A 



4J 



dk dk' G t (k,k')d(k-k' + Q) 
A, 



(27) 



which can as well be obtained from a self-consistent 
mean-field ansatz. For zero temperature, perfect nest- 
ing and a constant density of states, the gap behaves as 
A oc e~ 1 /( 4Jp °). Since we now have 



J xy =4JA/A . 



(28) 



J xy diverges for Ao , A — > and can thus be interpreted 
as a Goldstone vertex. J z , however, remains finite in the 
infrared and therefore corresponds to radial fluctuations 
of the staggered magnetization. The divergence of J xy for 
vanishing Ao also directly follows from the Ward identity 
Eq. [TJ which simplifies to 

A — A = — A J xy B-j-,4, . 



IV. GENERALIZED, ANISOTROPIC 
MEAN-FIELD MODEL 

In this section, we consider an interaction term of xxz- 
type 



S 



dk Vv (k) [ik - ek) ipr (k) + ^ S Q • S_ Q 



aJ 



S z n Sl 



(29) 



which breaks the SU(2)-symmetry for a non-vanishing 
anisotropy a. Eq. [29] interpolates between the isotropic 
case and an Ising-like interaction at a — ¥ oo for positive a 
or an interaction of xy-type for —1 < a < 0. The former 
will be referred to as the Ising regime and the latter as 
the xy regime in the following. The RG flow equations 
can be parametrized in exactly the same way as for the 
mean-field model in Eq. [T3] Again this leads to the flow 
eauations l22j|241 but now with the initial conditions J z = 
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FIG. 4. Flow of the gap A for a — 0.2, zero temperature, 
perfect nesting and a constant density of states po (W = 2/po, 
J = 0.09/po). 
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FIG. 5. Flow of J z for a = 0.2, zero temperature, perfect 
nesting and a constant density of states po (all other param- 
eters as in Fig. [4|. 



2J(1 + a), J xy — 4J and A = A in the ultraviolet. The 
formal solution of the flow equations now reads as 



J, = 



2J(l + a) 



1 + 47(1 + a) S t)t 
4J 



-4J(1 + a) %A. 



(30) 

(31) 
(32) 



Note that the Ward identity Eq. Q] has been derived for a 
spin-symmetric interaction and thus no longer holds. So 
a constraint that forces J xy to diverge in the limit Ao — > 
is now absent: for non-zero positive a, the Goldstone 
modes are gapped. As a consequence, the denominator 
of Eq. [31] no more needs be zero for Ao, A — > 0, but may 
now take on finite values. 

We now discuss the case of relatively weak anisotropics 
(a = ±0.2) for J > both in the Ising and in the xy- 



FIG. 6. Flow of J xy for a = 0.2, zero temperature, perfect 
nesting and a constant density of states po (all other param- 
eters as in Fig. [4]). The curves for different values of the bare 
gap Ao become discernable below the critical scale. 
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FIG. 7. Flow of A according to Eas.l22landl24lfor a = -0.2. 
zero temperature, perfect nesting and a constant density of 
states po (all other parameters as in Fig. |4]). Below X xy « 



10 / po the results are unphysical. 



regime. In Figs. |4]|6] the flow of the three couplings is 
depicted for a = +0.2. The flow of J z in Fig. [5] again is 
peaked at some critical scale A cr j t , which is now enhanced 
as can be seen from Eq.JJOl At this scale, the gap in Fig.0] 
starts to grow significantly and reaches a saturation value 
in the infrared. 

So far the results only quantitatively differ from the 
case of isotropic interaction. The crucial difference occurs 
in the flow of J xy . Inserting the modified gap equationl32l 



into Eq. [21] yields 



4J(1 



a 



A 



(33) 



which gives finite values in the limit A, A — > 0. 

We now consider the flow in the xy-regime with an 
anisotropy factor a = —0.2. Therefore the spins should 
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FIG. 8. Flow of J z according to Eqs.l22landl24lfor a = -0.2. 
zero temperature, perfect nesting and a constant density of 
states po (all other parameters as in Fig. Below X xy « 
8 ■ 10 _3 /po the results are unphysical. 



X z < X X y, where A starts to grow significantly. Note 
however that the formal solution of the flow equations [2^1 
for J z and [M] for A become unphysical at scales A < X xy , 
as the instability correspond to a staggered magnetiza- 
tion in the xy-plane. Our seed field, however, was sent to 
zero in the x- and y-direction before the thermodynamic 
limit was taken, which is not the physical order of the 
limits for spontaneous symmetry-breaking. Of course, 
with some additional effort, we could now allow for a 
xy-symmetry-breaking component of the self-energy in 
straight analogy to what we did so far for the z-direction. 
This would allow us to follow the flow further through the 
xy-instability. At the level presented here, we can infer 
that the formalism with symmetry breaking only in the 
^-direction is faithful enough to allow for additional in- 
stabilities. This is important in order to preserve the 
unbiased character of the RG approach. 



CHANNEL DECOMPOSITION 
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FIG. 9. Flow of J xy according to Eq. [33] for a = —0.2, zero 
temperature, perfect nesting and a constant density of states 
po (all other parameters as in Fig. [4j. At scales below the 
sign change of J xy the results are unphysical. 



be aligned in the xy-plane rather than in the z-direction 
in the ground state. Since our formalism does not al- 
low for a staggered in-plane magnetization, we should 
encounter an instability corresponding to spins ordered 
along a direction in the rry-plane. The flow of the gap 
and J z does not seem to differ much from the previous 
case. From Eq.[331 we however find that for — 1 < a < 0, 
J xy undergoes a sign change at SOIXLG SCcll© as shown 
in Fig. [9] Of course, this divergence forbids direct in- 
tegration of the flow equation [531 down to the infrared. 
Instead, a flow to strong coupling is observed. According 
to Fig.O this instability occurs at about X xy = 10 — 2 / po- 
At this scale, both A and J z remain finite. Formally, 
the flow of J z and A (Eqs. [2U and Y2M can be continued 
below X xy . The result obtained is similar to the previous 
cases (see Figs. [7] and [5]). Again J z has peak at a scale 



So far, we only have considered the RG flow for mean- 
field models. In realistic models showing an antifcr- 
romagnetic instability of itinerant electrons, however, 
transfer momenta other than the ordering vector Q play a 
role in the particle-hole channel, which, in addition, cou- 
ples with the particle-particle channel. If one studies the 
behavior of the effective interaction for wave vector trans- 
fers near Q, one will find that there is a sharp peak at Q. 
Similar statements hold for the frequency dependence, 
here a sharp peak forms around zero transfer frequency. 
As one will have to resort to a numerical evaluation of the 
flow equations and the corresponding wavevector and fre- 
quency summations in most cases, such sharp peaks pose 
a challenge to the usual discretization schemes. In the 
general case, these sharp structures can occur for both 
of the two possible wavevector/frequency transfers in the 
interaction t = k\ — k% and u — ki — ks, and for the 
total wavevector/frequency s = ki + &2- Hence the inter- 
action as function of k\ , ki and k% can possibly depend 
strongly on all these three variables. Here a vertex de- 
composition proposed by Husemann and Salmhofe r 14 i 15 
for the wave vector dependence and by Karrasch et aliSI 
for the frequency dependence can serve as a simplifying 
approximation. The idea is to write the interaction as a 
sum of (at least) three functions. One of these functions 
is used to capture the strong s-dependence of the inter- 
action, one is used for the i-dependence and on for the 
it-dependence, as will be detailed below. Then, instead 
of having one function that depends strongly on three 
variables, one has three functions that depend strongly 
only on one variable. Hence the discretization effort can 
be reduced from iV 3 down to 3iV. The idea is based on 
the observation that the second order corrections for con- 
stant bare interactions fulfill this decomposition exactly. 
One may hope that higher orders are still approximated 
well. 

So far, this decomposition has been proposed and 
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tested for the symmetric phase, which gave sufficient ev- 
idence for the validity of the underlying approximation. 
Recently, an extension to the superfluid phase has been 
given by Eberlein and Metzne r ' 19 . In this section, we 
propose a channel decomposition in the phase without 
spin-rotational symmetry, which should make a numeri- 
cal solution of the resulting flow equations tractable also 
in this case. 

For commensurate AFM, the renormalized interaction 
is only invariant under translations by an even number of 
sites. In momentum space, this means that the coupling 
functions are of the form 



V x (ki,k 2 ,k 3 ,k4,) = V x (ki,k 2 ,k 3 )6(ki + k 2 -l 
+ V x c (k u k 2 ,k 3 ) 5 Oi + k 2 -k 3 -k 4 + Q) 



with a momentum-conserving part V x an d a non- 
conserving part V x c , which is generated during the flow. 
In Nambu representation with spinors Eq. 1191 compo- 
nents W x " of the interaction with an even number of 
equal Nambu indices Si — ±1 correspond to V x while the 
components with an odd number of equal Nambu indices 
correspond to V x c . Thus we have 

W x (K 1 ,K 2 ,K 3 ,K i )=S {ki} 

V x (xi,X2,X 3 ) J2i s » even 
V x c {^,M 2l K 3 ) Odd ' 

with physical momenta Xi — ki + (1 — Si)Q/2 and 
Ki = (ki,Si). In this formula, the momenta are re- 
stricted to half the Brillouin zone, and therefore 6^} 
ensures momentum conservation up to multiples of 7r. We 
now decompose the three coupling functions from Sub- 



sec. |TTX]as follows. Renormalizations of equal-spin inter- 
actions Wf and W± can be regarded as a sum of triplet- 
superconductivity contributions and a spin-dependent 
CDW term 

W T {K 1 ,K 2 ,K 3 ,K 4 )=S {ki} \u^{k u k 2 ,k 3 ) 



{Si} 
tSC.T 

CDW. 

{Si} 
CDW,r 



(fci + fc 2 , ki, k 3 ) 

T (k 3 - k 2 ,ki,k 2 ) 
(ki - k 3 ,ki,k 2 ) 



where U T stems from the bare interaction and where 
s = (si, s 2 , S4, S3). The coupling function Wfi, in con- 
trast, is renormalized by a particle-particle coupling- 
function, which may contain triplet- as well as singlet- 
superconductivity terms, and magnetic contributions 
corresponding to S% + Sy and iS| terms in the interaction 

W^K^K^K^) =6 { h} [u\f{k u k 2 ,k 3 ) 

+ $ scti( fe i + k 2 ,k u k 3 ) 



SCf4 

*S, } (*i- 



- k 2 ,ki,k 2 ) 
k 3 ,k-L,k 2 



The right-hand sides of the flow equations Eqs. CTlUlnow 
read as follows. We use 

L%A{p,q) = d x [G%»(p) G s T ^(q)} 

as a short-hand notation for the loops and the mea- 
sure dip indicates that the respective momentum integral 
only runs over the folded Brillouin zone. In the particle- 
particle channel, one obtains 



S&VCU, q') = \ I'l' PW ' ['!■ I ~ q,P, I-P) W:'^ S3 ' s \l-p,p, q', I - q') (p, l-p), 
{<} 

Hc } n(hq,q') = -E fd' P W;i S2 ' s[ > s Hq,l- q,p,l-p)W;i' s *' S3 ' S4 (p,l-p,q',l- q') L^f (p, I - p) . 
{<} 



The particle-hole channel gives 



<i>cmv t ( J > ?> q') = - E f d 'p w t" S2 ' S3 ' s[ (i + p,q'J + q':P) w; 1 ' s ^ Si ( q , P ,i + P ,q-i) l& \ p ,i + P ) 
{<} 

- fd'pW;i' s '^(q', l+pj + q',p) W;i s '^' s Hq,P, q-Ul + P) 4j } (p, l + p), 
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|0> 9, </) = - J2 [d'pWf' S2 ' S3 ' s i(l + P, q', l + q',p) W^ s '^ Si (q,p, l+p,q-l) L}f(p, l+p) 
{<} 

- I d 'P W U S2 ' S[ ' S3 {l+P,q',P,l + q') W;i' SUS '*' Si (p,qj+p,q-l) z\f (p, l + p), 



= - E [d'pW;'l' S2 ' S3 ^(l+p,q',l + q',p) W;i' s '*' s *> Si (q,p,l+p,q- l)L$(p,l+p) , 
{<} 



{<} 

-J2 Up W^' S2 (l+p,q',l + q',p) ' S '*< S3 ' S '° (q,p,q-l,l+p) lJ* S } (p, / + p) . 
{*«} 



This channel decomposition should allow for efficient fur- 
ther parameterization giving rise to numerically tractable 
flow equations. Nevertheless, implementing and testing 
this scheme represents a longer effort, and at this point 
we cannot present any specific results. 



VI. SUMMARY AND OUTLOOK 

In this work we have parameterized the effective action 
and given fRG flow-equations for spin- 1/2 fermions in 
phases with broken spin-SU(2) invariance. One then has 
to deal with three coupling functions for the interactions 
instead of one in the symmetric phase. The flow of the 
AFM mean-field model can be captured by three running 
couplings corresponding to radial and Goldstone compo- 
nents of the interaction and the gap. As one may have ex- 
pected, these couplings behave as in previous studies^^ 
of other mean-field models. The radial coupling J z shows 
a pronounced hump at the critical scale, where the gap A 
starts to grow significantly. While the radial coupling and 
the gap saturate to finite infrared values for small seed 
fields, the Goldstone coupling J xy diverges, reflecting the 
masslessness of the Goldstone modes. This can as well be 
seen from the global SU(2) Ward-identity. The RG flow 
equations reproduce the RPA result for the interaction 
and the self-consistent gap equation exactly. This sug- 
gests that we have indeed found a sensible parametriza- 
tion. A generalization to a symmetry-breaking interac- 
tion of xxz-type is straightforward. In the Ising regime, 
the Goldstone modes are then gapped, which renders all 
couplings finite in the infrared. For negative anisotropy 
parameters of the bare interaction, i.e. in the xj/-regime, 
J xy flows to strong coupling. This instability corresponds 
to the ordering of spins in the xy-plane. However, this 
type of spontaneous symmetry breaking is not included 
in the formalism presented here (but could be included 



I 

with some additional work). This example shows that 
the built-in symmetry-breaking in the ^-direction going 
along with a gapping of the Fermi surface does not pre- 
vent the physically reasonable instability toward ordering 
in the xy-plane. Hence, at least on this level, triggering 
symmetry breaking in a specific direction does not re- 
move the usual RG advantage of providing an unbiased 
description. 

For an isotropic bare interaction, this formalism can 
also be used for proceeding beyond the level of a reduced 
mean-field model. Then the usual advantage of RG ap- 
proaches, which is to include all the different one-loop di- 
agrams, comes into play and can lead to significant devia- 
tions from the mean-field picture. However, on the tech- 
nical level various complications arise: Particle-particle 
diagrams lead to a non-trivial momentum and frequency 
dependence of the coupling functions. Since this has to 
be resolved numerically or by a suitable ansatz, we expect 
that this drastically increases the computational cost. 
Hence a channel decompositio n 14 ! 15 appears to be in or- 
der. This decomposition replaces an interaction that may 
depend strongly on three wave vectors and frequencies by 
(at least) three functions that depend strongly only on 
one wave vector sum or transfer. This reduces the nu- 
meral effort considerably. For commensurate AFM, the 
symmetry breaking one-particle term violates momen- 
tum conservation and non-conserving contributions to 
the interaction are generated in the flow. In the commen- 
surate case, the Nambu formalism helps to facilitate our 
notation. We therefore have derived channel-decomposed 
flow equations in Nambu representation. The next step 
will indeed be the implementation of this scheme for non- 
mean-ficld models. As this is still a larger numerical 
task, we refrain from embarking on this here and only 
describe the theoretical foundation. In solving the flow 
equations for the more general case, one may proceed as 
for a fermionic superfluid 1 ^. However, especially in the 
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case of AFM, an additional feature may be of interest. In 
Ref. [l9l the momentum dependence of both the fermion- 
boson vertices and the gap was assumed to correspond to 
a simple angle-independent s-wave, which is a reasonable 
approximation for the attractive Hubbard model. (How- 
ever, the approach pursued in that work can in principle 
account for more complicated boson-fermion vertices.) 
In the repulsive Hubbard model at small but non-zero 
second-neighbor hopping the situation is different in two 
ways. For one thing, the particle-particle contributions to 
the interaction may then develop some structure that has 
an impact on the fermion-boson vertices in the particle- 
hole channels, e.g. on their magnitude and dependence 



on the fermionic wavevectors. It may therefore be ad- 
vantageous to track the deformation of fermionic form 
factors in the flow. Furthermore, a reconstructed Fermi 
surface might still give rise to a Cooper instability inside 
the AFM symmetry-broken phase, which would account 
for a homogenous coexistence of non-zero AFM and su- 
perconductivity order-parameters. In principle, this ef- 
fect should be captured in our formalism, but its demon- 
stration remains an issue of future research. 

We thank A. Eberlein, A. A. Katanin, W. Metzner, 
and M. Salmhofer for useful discussions. This project 
was supported by the German Resarch Foundation DFG 
via FOR 723. 
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